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Neutrinos of astrophysical origin are messengers produced in stars, in explosive phenomena like 
core-collapse supernovae, in the accretion disks around black holes, or in the Earth's atmosphere. 
Their fluxes and spectra encode information on the environments that produce them, but they 
are modified when neutrinos traverse a medium. Here we summarize our current understanding of 
neutrino flavour conversion in media and emphasize the importance of this domain for astrophysical 
observations. We discuss examples of the fundamental properties that astrophysical neutrinos have 
uncovered, or might reveal in the future. 



I. INTRODUCTION 

R. Davis' pioneering measurement of electron neutri- 
nos emitted by the sun has opened the era of neutrino 
astronomy [1 1. A significant deficit, compared to the 
standard solar model predictions, was rapidly measured, 
triggering several decades of experiments to investigate, 
if unknown neutrino properties, or solar model uncertain- 
ties, were at its origin. Electron antineutrinos produced 
by a massive star were first detected during the SN1987A 
explosion. A burst was first seen in Kamiokande, the elec- 
tron direction pointing to the Large Magellanic Cloud at 
50 kpc from the Milky Way. The occurrence probabil- 
ity of having 9 events per 10 seconds was determined to 
be less than 5.7 10 -8 [2]. Altogether the Kamiokande, 
1MB, Baksan and Mont Blanc detectors observed about 
twenty events [3HS], in the 40 MeV energy range during 
13 seconds (Figure 1). Their angular dependence and en- 
ergy distribution are reasonably consistent with expecta- 
tions from core-collapse supernova simulations (see e.g. 
[7] for a recent analysis). The detection of these events 
constitute the first experimental confirmation of super- 
nova simulation predictions. In 2002, R. Davis and M. 
Koshiba were awarded the Nobel Prize, for "pioneering 
contributions to astrophysics, in particular for the detec- 
tion of cosmic neutrinos", shared with R. Giacconi for 
"for pioneering contributions to astrophysics, which have 
led to the discovery of cosmic X-ray sources" . The puzzle 
of the "solar neutrino deficit" was solved by the Super- 
Kamiokande experiment with the discovery of neutrino 
oscillations in 1998 8 . A deficit of upgoing athmospheric 
muon neutrinos traversing the Earth was observed, com- 
pared to downgoing ones, that was shown to be consis- 
tent with the hypothesis that upgoing muon neutrinos 
convert into tau neutrinos. This discovery has fundamen- 



tal implications for high-energy physics, astrophysics and 
cosmology. For example, neutrinos from the sun give us 
direct confirmation of the energy production mechanisms 
in stars (see the seminal works [91 [TU]). 

The neutrino oscillation phenomenon occurs if neutri- 
nos have non-zero masses and mixings. The possibil- 
ity to have v-v oscillations was first pointed out by B. 
Pontecorvo in analogy with Kq — Ko mesons [TT] . Neu- 
trino oscillations require that the flavor v a and the mass 
Vi basis are related by v a = U*.^ with U being the 
Maki-Nakagawa-Sakata-Pontecorvo (MSNP) unitary ma- 
trix [T2]. In the three-flavor framework U = VD with a 
possible parametrization given by (T3] : 
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where Cy = cos Oij and Sjj = sin 9^ and D = 
diag(e~ l< ^ 1 , 1, e~ l< ^ 2 ). The V matrix is analogous to the 
Cabibbo-Kobayashi-Maskawa matrix, in the quark sec- 
tor, although the mixing angles are there measured to be 
small. The <pi , </> 2 are two extra phases that appear in the 
Majorana case. Therefore the MNSP matrix depends on 
three mixing angles, two Majorana-type and one Dirac- 
type phases. (Here we will not discuss the Majorana 
phases, since they do not influence neutrino oscillations 
in vacuum and in matter.) If the Dirac-type phase is 
non-zero, the MNSP matrix is complex, introducing a 
difference between the neutrino oscillation probability, 
and the corresponding ones for antineutrinos, implying 
CP violation in the lepton sector (for a neutrino physics' 
overview, see the recent comprehensive book 14J). 

Oscillations in vacuum is an interference phenomenon 
among the matter eigenstates, which is sensitive to 
their mass-squared differences and to the mixing an- 
gles. In the two-flavour framework, the oscillation ap- 
pearance probability for relativistic neutrinos is given 



2 



by P(y a -> i/ja) = sin 2 (26l)sin 2 (Amf 2 L/4£;), E being 
the neutrino energy, L the source-detector distance and 
Am 2 ,! = mj — m\. While the oscillation amplitude de- 
pends on the mixing angle, the squared-mass differences 
determine the oscillation frequencjQ In the last decade, 
reactor, accelerator and solar experiments have precisely 
measured the two mass-squared differences to be Am^ = 
(7.50 ± 0.20) 10" 5 eV 2 , |Am§ 2 | = (2.32 + 0.12) 10" 3 eV 2 , 
as well as the mixing angles sin 2 (26*12) = 0.857 ± 0.024, 
sin 2 (26>i 3 ) = 0.098±0.013 Q3]. The third angle has been 
recently measured by the T2K [15], the Double-Chooz 
[US] , the RENO [T7J and the Daya-Bay collaborations 
giving sin 2 (26> 13 ) = 0. 092 ± 0. 016(stat) ±0. 005 (syst) [IS]. 

While the experimental progress is impressive, numer- 
ous features remain unrevealed. The mechanism that 
gives to the neutrino a mass is unknown. Depending on 
the neutrino nature, it might require the introduction of 
a right-handed neutrino singlet (that does not couple to 
the gauge bosons) or more complex mass mechanisms of 
the see-saw type(s), for example, which would ask for 
significant Standard Model extensions. There is, what is 
usually referred to, as the mass hierarchy problem. Since 
one of the Am 2 signs is unknown, there are two ways to 
order the mass eigenstates : Am 2 3 > (normal hierar- 
chy), or Am 2 3 < (inverted hierarchy). Information on 
the absolute mass scale is extracted from tritium beta- 
decay experiments. The current upper limit of about 2 
eV for the electron neutrino effective mass [13] will be 
significantly improved by the KATRIN experiment that 
is expected to reach 0.35 eV sensitivity [19]. Moreover, 
key open issues are the existence of leptonic CP viola- 
tion, of sterile neutrinos and the Dirac versus Majorana 
v nature. 

Interestingly, the ensemble of experimental data 
present a few anomalies that cannot be cast in the three- 
active neutrino framework. The neutrino flux measure- 
ment from intense static 37 Ar and 51 Cr sources in the 
GALLEX and SAGE experiments present an anomalous 
deficit of electron neutrinos (the "Gallium anomaly"). 
A recent statistical analysis finds a significance at 3<r 
[20] . After the controversial LSND results [21], Mini- 
B00NE anti-neutrino and neutrino oscillation data com- 
bined show an excess of events at 3.8(7 [35] [53]. Besides 
a recent re-evaluation of the electron anti-neutrino flux 
from reactors has shown a shift in the flux renormaliza- 
tion by 3 % compared to previous predictions. The re- 
analysis of all reactor experiments using this new flux has 
shown a a 3cr inconsistency with the standard oscillation 
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1 For three flavors only two independent Am 2 can be built. To 
accommodate an extra Am 2 a fourth mass eigenstate is required. 



FIG. 1. Electron anti-neutrino events measured by the 
Kamiokande, Baksan, 1MB and LSD experiments during the 
supernova 1987A, in the Large Magellanic Cloud [2]. 



framework (the "reactor anomaly") 24 . These unex- 
plained features might point to new physics and require 
one sterile neutrino, non-standard interactions or CPT 
violation. At present none of the proposed explanations 
is capable of providing a comprehensive understanding. 

Future experiments using man-controlled sources will 
address unknown neutrino properties and the possible 
explanations of the anomalies. Exploring the Dirac CP 
violating phase requires long-term accelerator measure- 
ments employing either established or novel techniques, 
such as super-beams, beta-beams [25] , neutrino factories 
(see e.g. [Ml [37]), or projects like the decay-at-rest based 
DAEDALUS [28]. With upgrades of T2K and NOz/A ex- 
periments a small fraction of the 6 values can also be 
investigated at 3a level [29]. Majorana- type phases and 
the neutrino nature can be determined by searches for the 
lepton-number violating neutrino-less double-beta decay. 
Earth matter effects in long-baseline experiments can be 
exploited to identify the hierarchy, as in e.g. [27]. Numer- 
ous projects are being conceived to test the hypothesis 
of the existence of sterile neutrinos based e.g. on the use 
of intense radioactive sources inside spherical detectors 
such as 144 Ce in Ref.^U] or 8 Li in Ref.[3T] (see [32] for a 
review) . 

Astrophysical and cosmological observations offer com- 
plementary strategies in these fundamental searches. Al- 
though indirectly, cosmological data will reach an un- 
preceded sensitivity on the sum of the neutrino masses 
[33] [34] and maybe be sensitive to the hierarchy. Athmo- 
spheric neutrinos in experiments like PINGU [35 in Ice- 
Cube, or ORCA [35] in ANTARES, might offer a unique 
tool to determine the mass hierarchy. The neutrino mass 
ordering also leaves characteristic imprints in the time 
and energy signals of core-collapse supernova neutrinos, 
if an (extra)galactic explosion occurs. This would be seen 
by a network of detectors, such as Super-Kamiokande, 
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KamLAND, Borexino, or by dedicated supernova obser- 
vatories like HALO and LVD; while some constitute the 
"SNEWS: SuperNova Early Warning System" aiming at 
alerting observers if a galactic supernova explosion occurs 
[3"T] . About a thousand events can be collected if a super- 
nova explodes in our galaxy. While such phenomena are 
rare, the integral over cosmological redshift of supernova 
neutrinos produces the yet unobserved diffuse supernova 
neutrino background. The sensitivity for its discovery 
could be reached with improved technologies [35] or with 
large-size detector^ like MEMPHYS or Hyper-K (440 
or 770 kton water Cherenkov) and GLACIER (100 kton 
liquid Argon) [39J that would collect about 350 and 60 
events over 10 years, respectively |40) . Key information 
from the diffuse supernova background detection could 
be extracted on the supernova dynamics, on the star for- 
mation rate and on unknown neutrino properties (for a 
review see [4T1 142]). 




II. NEUTRINO FLAVOUR CONVERSION IN 
MEDIA : STATUS AND OPEN QUESTIONS 



FIG. 2. The MSW effect : Schematic behavior of the neutrino 
matter eigenvalues (solid lines) as a function of the matter 
number density of a medium. The eigenvalues are pushed 
far apart because of the presence of mixings, at the MSW 
resonance location. (The dashed lines shows the eigenvalues 
in absence of mixings.) In the case of an adiabatic evolution 
the matter eigenstate vb does not mix with the va at the 
resonance and stays on the higher branch up to the surface of 
the star . An electron neutrino, born as vb emerge as vb, for 
an adiabatic evolution, which produces an electron neutrino 
deficit in a i/ e sensitive detector on Earth. 



It is now established that the neutrino solar deficit is 
due to the Mikheev-Smirnov-Wolfenstein (MSW) effect 
[43} I44j . a resonant conversion phenomenon occurring 
when neutrinos interact with the matter composing a 
medium. In two-flavourSjthe MSW resonance location is 
identified by the relation^] V2G F n e = (Am 2 cos 20)/ 2E, 
with Gf the Fermi coupling constant, n e the electron 
number density |14j . At such a location, a conversion 
from the electron to the muon (and tau) flavours can take 
place, whose efficiency (or adiabaticity) , also depends on 
the mixing angle and mass-squared difference values on 
one hand, and on the matter number density profile on 
the other (see Figure 2). According to the squared-mass 
difference sign, the MSW effect occurs either in the elec- 
tron neutrino channel, or in the electron anti-neutrino 
one. Since R. Davis' experiment, numerous solar exper- 
iments mainly sensitive to the electron flavour have pre- 
cisely measured the solar neutrino flux. Using charged- 
and neutral-current neutrino scattering on heavy water, 
the SNO experiment has showed that the measurement of 
the total neutrino flux is consistent with the predictions 
of the standard solar model : solar electron neutrinos 
convert into the other active flavours [IS]. Moreover the 



2 Such detectors are about 20 times larger than Super- 
Kamiokande. Their goal is cover an interdisciplinary program 
including the detection of supernova neutrinos, leptonic CP vio- 
lation and proton decay searches. 

3 For antineutrinos the r.h.s. of the relation has a minus sign. 



reactor experiment KamLAND has identified the Large 
Mixing Angle (LMA) solution, by observing reactor elec- 
tron anti-neutrino disappearance at an average distance 
of 200 km, definitely determining that neutrino flavour 
conversion in matter produces the solar deficit |46j . Fig- 
ure 3 shows the comparison of solar experiments' results 
with the LMA prediction. Nowadays the MSW effect 
constitutes a reference mechanism in the study of neu- 
trino flavor conversion in media. 

Another interesting case is represented by neutrinos 
traversing the Earth. Its matter density profile undergoes 
a rapid change from the mantle to the core at higher den- 
sity [H]. It was first pointed out in [32] and in [SOj that 
atmospheric neutrinos traversing the Earth would change 
in flavour. Besides the MSW effect, another neutrino 
flavour conversion mechanism can occur, named "para- 
metric resonance" [5T] or " neutrino oscillation length res- 
onance" [52] . This effect is an interference effect due to 
the mantle-core-mantle change in the Earth matter den- 
sity profile, that can enhance the oscillation probabilities. 
For example, depending on the trajectory (or azimuthal 
angle) and energy, atmospheric neutrinos in the few GeV 
energy range can have an MSW resonant conversion in 
the core, or in the mantle and in the core, and for the 
core-crossing trajectories also feel the parametric reso- 
nance. The MSW and parametric resonance effects have 
been investigated e.g. in [53] in the subdominant atmo- 
spheric v e — > in Super-Kamiokande. A vaste literature 
concerns matter effects in atmospheric neutrinos (see e.g. 
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FIG. 3. Solar neutrinos : Electron neutrino survival prob- 
ability as a function of the neutrino energy. The results cor- 
respond to pp, pep, 7 Be, S B solar neutrinos. The gray band 
shows the predictions assuming that the MSW effect occurs 
in the sun with the LMA solution (the band corresponding to 
the la range of the mixings parameters) [47] . 



[54"] and references therein). 

A variety of novel flavour conversion mechanisms has 
been identified in stars with more than 6-8 solar masses 
- the O-Ne-Mg and iron core-collapse supernovae. Their 
evolution ends with an explosion where 99 % of the grav- 
itational energy is released as neutrinos of all flavours 
in the 100 MeV energy range in a ten seconds' burst, 
leaving either a neutron star, or a black hole. Since 
the matter number density of these stars is very large, 
the MSW effect occurs at two different locations, usually 
referred to as the High-resonance (H), associated with 
(#13, Am! 3 ), and the Low-resonance (L), corresponding 
to the (0 12 ,Am? 2 ) (see Figure 4) [55]. A 40 MeV neu- 
trino sees the H-resonance at a density of about 10 3 
g/cw 3 . The identification of the solar LMA solution tells 
us the value and the sign of Am 2 2 . The three neutrino 
mixing angles are also now precisely determined. Ac- 
cording to typical supernova matter density profiles, the 
neutrino passage in the L-resonance region is adiabatic, 
while the neutrino behaviour at the H-resonance encodes 
interesting information on the neutrino mass hierarchy. 
Therefore matter effects in supernova neutrinos can be 
exploited to learn about the hierarchy using signals from 
future explosions. This line of research is followed by nu- 
merous works proposing signatures that are specific to 
the mass hierarchy being normal or inverted. However, 
theoretical investigations of the last ten years have shown 
that the situation is more complex than the framework 
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FIG. 4. MSW effect in supernovae : Level crossing dia- 
gram for the neutrino mass eigenstates in a supernova for the 
case of normal mass hierarchy. The figure shows the matter 
eigenvalues evolution associated with the matter eigenstates 
v\ m , V2m, vz m , as a function of the electron number density n e . 
The two crossings correspond to the high-resonance (H) and 
the low-resonance (L) associated with the mixing parameters 
(#13, Amy) and (#12, Amf 2 ), respectively [55] . 

shown in Figure 3 of Ref . [5 5) , because one is facing an ex- 
plosive phenomenon with shock wave(s) and turbulence 
on one hand, and where 10 58 MeV are being released as 
neutrinos on the other. 

Ref. [56] pointed out that the shock wave passage in 
the MSW resonance region would leave an imprint on 
the time signal of the positrons emitted in inverse beta- 
decay, i.e. v e + p —> n + e + , the main detection chan- 
nel in Cherenkov and scintillator detectors. Ref. [57J 
has shown the presence of not only a front but also a 
reverse shock. The presence of shock waves engenders 
two important effects. First it makes the H-resonance 
temporarily non-adiabatic because of the steepness of 
the shock wave fronts (the adiabaticity depends on the 
derivative of the density profile). Second, multiple H- 
resonances appear since a neutrino of a given energy can 
meet the resonance condition several times. For exam- 
ple the electron neutrino survival probability in pres- 
ence of two resonances reads P ee — cos 2 (xi — X2) — 
sin 2%i sin 2\2 sin 2 ( Am 2 / AE)dx with xi an d X2 the 
matter angles mixing the two matter eigenstates at the 
resonance locations x\ and a; 2 , and Am 2 the matter 
squared-mass difference. The interference term oscillates 
with the neutrino energy and the resonance locations. 
Since the shock wave is moving, it will change the reso- 
nance locations, modifying the phase in such a term. This 
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produces rapid oscillations in the survival probability for 
a neutrino of a given energy. The conditions for such 
oscillations termed "phase effects" are that the flavour 
evolution is semi-adiabatic at the resonances and that 
the matter eigenstates stay coherent at the resonances 

ED- 

Turbulence is another characteristic of these explosive 
environments. Its effects are studied e.g. in [59"rl61| . 
The presence of matter density fluctuations can intro- 
duce numerous locations where the resonance density is 
met. The effect on flavor conversion is therefore of the 
same kind as for the case of shock waves in the sense that 
turbulence produces phase effects, averaging the neutrino 
oscillation probability between the adiabatic and non- 
adiabatic solution [62] (see [63] for a review). Extensive 
calculations with matter density fluctuations extracted 
from multi-dimensional supernova simulations remain to 
be done. Clearly, to fully capture several flavour con- 
version features in media, that are interference effects, 
one has to evolve the neutrino amplitudes and not the 
neutrino probabilities. The latter procedure was often 
followed in the past, e. g. in predictions with the factor- 
ization hypothesis where one factorizes the probabilities 
for flavor conversion at the H- and L-resonances. 

The implementation of the neutral-current neutrino- 
neutrino interaction is another important progress in sim- 
ulations of flavour conversions in supernovae. Such con- 
tributions are small in the sun, whereas the large neutrino 
number density renders them important in such enviro- 
ments. Pantaleone first pointed out that the neutrino 
self-interaction would introduce a non-linear refractive 
index |64| while the first numerical simulations showed 
the appearance of new phenomena [55]. Ref.[B"6"] pointed 
out a significant impact on the r-process nucleosynthesis 
and Ref. [57] showed the important effects on the neu- 
trino fluxes with phenomenological implications, trigger- 
ing intensive investigations (see |68] for a review). Pre- 
dictions are obtained by solving a system of coupled stiff 
and non-linear equations for the neutrino density matrix 
(or amplitudes or polarization vectors). Three conver- 
sion regimes are indentified. At first the neutrino flavor 
conversion is frozen (the "synchronization"), then bipo- 
lar oscillations are triggered by the presence of flavour 
instabilities. Such a regime can be described as a "fla- 
vor" [69] or a "gyroscopic" pendulum [70]. Then, full 
or no flavor conversion occurs depending on the neu- 
trino energy, whose underlying mechanism, in the comov- 
ing frame, has an MSW-like behaviour [71], which can 
be interpreted as a "magnetic resonance" phenomenon 
[72) . This mechanism produces sharp changes of the neu- 
trino fluxes, called "spectral splits" : the electron and 
muon (or tau) neutrino fluxes swap above a critical en- 
ergy (the "split" energy). Note that these effects are 



also referred to as "collective effects" because many fea- 
tures can be captured by following the neutrino ensemble. 
These mechanisms explain the flavour modifications ap- 
pearing in simulations, while the picture becomes more 
complex depending on the neutrino luminosity ratios at 
the neutrino-sphere]^] on the neutrino properties and on 
the detailed implementation of the geometry. For ex- 
ample, collective effects can be modified by the presence 
of large matter densities that can introduce decoherence 
[73"] . Calculations based on realistic matter density pro- 
files from one-dimensional supernova simulations show 
that they might be suppressed during e.g. the accre- 
tion phase of the supernova explosion |74) . In this case 
the situation becomes simpler since the flavour evolu- 
tion picture reduces to the MSW effect, but with the 
modifications introduced by the shock wave effects and 
turbulence. 

Impressive progress has been achieved in unravelling 
mechanisms and conditions for flavour conversion of 
core-collapse supernova neutrinos. Still, serious work is 
needed to come to a definite and comprehensive under- 
standing and to establish the impact on observations. 
Up to now, the neutrino transport in the dense super- 
nova region is implemented in detail in supernova sim- 
ulations that typically account for a good angular, or 
energy neutrino distribution (but not both), and do not 
include mixings |75j . The neutrino energy spectra and 
fluxes at the neutrinosphere from such simulations are 
taken as initial conditions in the flavour evolution stud- 
ies, that evolve v up to the star surface, to predict signals 
in observatories on Earth. The corresponding evolution 
equations are based on the mean-field approximation : 
one neutrino (or anti-neutrino) is evolved at a time in a 
background of matter and of (anti)ncutrinos, that acts 
on the "test" particle through the mean-field created by 
them. It is still an open question to know, if and un- 
der which conditions, the present treatments describe 
the transition region that is Boltzmann treated to the 
one that is mean-field described, in a appropriate way; 
or if one needs more realistic geometrical descriptions, 
or extended equations implementing many-body correc- 
tions. For example, Ref . [76] has pointed out the need 
for a more realistic geometrical description of the neu- 
trino emission at the neutrino-sphere (Figure 5). Ref. [77] 
has derived corrections beyond the mean-field approx- 
imation using a coherent-state path integral approach. 
Ref. [75] has used the general many-body framework of- 
fered by the Bogoliubov-Born-Green-Kirkwood-Yvon hi- 
erarchy to derive both the mean-field and extended mean- 



4 The neutrinosphere is the region in the supernova where neutri- 
nos start free streaming (Figure 5). 



FIG. 5. Cartoon picture of neutrinos emitted at the neu- 
trinophere of a supernova, taken to be a sharp sphere, and 
interacting with each other |76] . 



field equations. The latter also introduce for the first 
time an abnormal mean-field corresponding to neutrino- 
antineutrino pairing correlations (Figure 6). So far ne- 
glected, these two-body correlations might impact flavour 
conversion and supernova observations. 

The theoretical developments described above also ap- 
ply to the study of the low energy neutrinos from the 
accretion disks surrounding black holes (AD-BH). In fact 
the conditions in the disks are very similar to the core- 
collapse supernova ones, so that such sites also offer con- 
ditions for example for r-process and vp element nucle- 
osynthesis (see e.g. jTHUSH])- Refs. [6"6"II5T] has studied e.g. 
the effect of the neutrino-neutrino interaction in these 
sites. Finally, neutrino flavour conversion phenomena 
are also important in the Early Universe, in particular at 
the epoch of Big-bang nucleosynthesis. Their description 
is based on the resolution of a Boltzmann equation for 
particles with mixings since collisions need to be imple- 
mented, as well as neutrino mixings, the coupling to the 
relativistic plasma and neutrino self-interaction whose 
importance was indeed first pointed out in this context. 
As for the r-process, interactions of electron neutrinos 
and anti- neutrinos with protons and neutrons, but also 
neutron decay, determine the proton-to-neutron ratio, a 
key parameter for primordial nucleosynthesis. The build 
up of the light elements starts at about 1 s after the Big- 
bang when the plasma temperature is around 1 MeV and 
is very sensitive to neutrino properties (see [HH IEH] for a 
review). For example, current cosmological observations 
from primordial nucleosynthesis and the matter power 
spectrum are compatible with about one sterile neutrino 
(see e.g. [El]). The upcoming results from the PLANCK 
experiment should tell us precisely how many extra de- 
grees of freedom cosmology allows. As for the direct 
observation of the cosmological neutrino background, it 
remains a challenge. Since such v are very cold, novel 
detection approaches are necessary. Attempts are still 



FIG. 6. The left figure represents neutrino-antineutrino pair- 
ing correlations, one of the interaction terms present in the ex- 
tended mean- field neutrino evolution equations in supernovae. 
The right figure shows the corresponding pairing mean-field, 
obtained by summing over the background particles states. 
Similar figures can be drawn for the usual electron and neu- 
trino background contributions to the normal mean-field [78] . 



off by several order of magnitudes from the required sen- 
sitivity. The one which appears the "closest" is the de- 
tection of cosmological neutrinos through their capture 
of radioactive nuclei, a process without threshold. First 
proposed by Weinberg |85) , this idea is applied in realistic 
extensive calculations in [86] (see also Ref. [87] ) . 



III. OPEN ISSUES AND ASTROPHYSICAL v 

One of the crucial open issues is the possible breaking 
of CP invariance. The recent measurement of a non-zero 
third neutrino mixing angle by the T2K [T5] , the Double- 
Chooz PI], the RENO [H] and the Daya-Bay collabora- 
tions [18] brings a much awaited result since CP violation 
from the Dirac phase affects oscillations only in the case 
of three families. Observing leptonic CP violation and 
neutrino-less double-beta decay would give important 
clues for the understanding of matter-antimatter asym- 
metry in the universe (see e.g. [EE])- An interesting ques- 
tion to explore is the possibility to have indirect manifes- 
tations of a non-zero Dirac-type phase in astrophysical 
environments. The authors of Ref. [89] have first investi- 
gated this question in the solar neutrino context showing 
that there are no CP violating effects in the v e survival 
probability in matter to the leading order in electroweak 
interaction while those from next-to-leading order are es- 
timated to be extremely small. Ref. [5D] has first pointed 
out possible effects in core-collapse supernovae coming to 
the conclusion that there are no CP violation effects in 
such environments. Ref.[9Ij has come to a result at vari- 
ance with Ref. 90 demonstrating that there can be CP 
violating effects in supernovae. By using general argu- 
ments, the factorization condition H(S) = H(8 = 0)S 



7 



is obtainecj^] with S = diag(l,l f e iS ). Such a relation 
gives a procedure to identify under which conditions such 
effects can arise : contributions to the neutrino Hamil- 
tonian that break such a condition engender CP violat- 
ing effects. For example, these arise because of radiative 
corrections in the Standard Model, or of non-standard 
interactions like flavor-changing interactions, that differ- 
entiate muon from tau neutrino interactions with mat- 
ter. This finding has been subsequently independently 
confirmed by Ref. [52] • A quantitative evaluation of the 
neutrino flux modification shows variations at the level 
of a few percent, both on a possible signal in an observa- 
tory on Earth, or on the neutron-to-proton ratio relevant 
for r-process nucleosynthesis [91] . An improved numer- 
ical simulation comprising neutrino coupling to matter 
and to neutrinos has shown that the non-linearity of the 
equations amplifies the CP violating effects from radia- 
tive corrections, for example, by several order of magni- 
tudes [93 . While flux modifications are still at the level 
of about 5-10 %, improvements of the simulations might 
introduce further amplification of such effects. A similar 
study has been performed in |94j in the Early universe 
context, where a non-zero Dirac-type phase is consid- 
ered for the first time. Within three-active neutrinos it 
is shown that the phase modifies the primordial Hclium-4 
abundance at most by 1%, which is within current sys- 
tematic uncertainties. 

The detection of neutrinos from a future (ex- 
tra) galactic explosion constitutes one of the strategies to 
identify the mass hierarchy. One possibility is to look for 
the presence of dips in the time signal due to the shock 
waves as in [57]. Ref. [55] also considers the uncertainties 
on the neutrino fluxes from supernova simulations. Note 
that these works do not implement the v — v interaction 
yet. In Ref. [96] the positron time signal associated with 
inverse beta-decay, in water Cherenkov and scintillator 
detectors, is predicted. The calculation treats neutrino 
self-interactions and shock-wave effects in a consistent 
way, using realistic supernova density profiles and prop- 
agating neutrino amplitudes (phase effects are properly 
taken into account). The result is shown in Figure 7. At 
early times, when the shock wave has not yet reached 
the MSW resonance region, the neutrino conversion is 
expected to be adiabatic. As the shock wave reaches this 
region, the flux becomes the non-adiabatic one, produc- 
ing either a dip or a bump, depending on the neutrino 
energy. For electron anti-neutrinos, the resonance con- 
dition is met for inverted hierarchy. If the hierarchy is 
normal, no resonant conversion occurs and the positron 
time signal presents an exponential behaviour. Such a 



5 H being in the so-called T23 basis 




t[s] 



FIG. 7. Positron time signal associated with inverse beta- 
decay for a galactic explosion at 10 kpc per unit tonne in a 
detector. The different curves are for 10 (solid), 15 (dashed), 
19 (dash-dotted) and 29 (dot-dot-dashed line) MeV positron 
energies. The dips or bumps correspond to the passage of the 
shock wave in the MSW resonance region. They are present 
in the case of inverted mass hierarchy |96| . 

signature could already be seen in Super-Kamiokande if 
a galactic supernova happens tomorrow, and be distin- 
guished by the exponential at 3.5er (lc) for the bump 
(dip). Note that this signature holds in the electron neu- 
trino channel, if the mass hierarchy is normal and could 
be seen in a liquid argon detector for example, mainly 
sensitive to v e , or through scattering on oxygen in water 
Cherenkov and carbon in scintillator detectors. Various 
other signatures of the mass hierarchy have been inves- 
tigated in the literature (see e.g. [57]). In Ref. [55] the 
signal of the first 100 ms (early post-bounce) of the explo- 
sion is investigated in IceCube. This has the advantage 
of being simpler than the one of the later explosion stages 
since it is not affected by the neutrino-self-interaction, or 
the shock- waves (and turbulence). 

Complementary information would be obtained by the 
observation of (extra) galactic supernova explosion(s) and 
the detection of the diffuse supernova neutrino back- 
ground. Clearly, the number of degeneracies reduces by 
the combination of data from detectors based on different 
technologies and with various energy thresholds for the 
neutrino and anti-neutrino detection channels. Figure 
8 shows an example by considering the one- and two- 
neutron detection channels in a lead-based detector like 
HALO [55J. Note that a software - SNOwGLoBES - is 
now available to compute interaction rates for supernova 
neutrinos in common detector materials (100[ 1101] . 
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FIG. 8. One- and two-neutron emission event rates for a 
supernova at 10 kpc in the HALO phase-II detector (1 kton 
of 208 Pb). The different curves take into account the uncer- 
tainties in the supernova muon (tau) neutrino fluxes at the 
neutrino-sphere as well as the possible hierarchy (#13 is now 
measured) [99] . 



As for observations using the v e detection chan- 
nel, they would benefit from precise measurements of 
neutrino-nucleus scattering. Only for deuteron the cross 
section is known at the level of few percent, while for 
heavier nuclei uncertainties are typically at the level of 
few tens of percent. In the future facilities producing neu- 
trinos in the 100 MeV energy range, such as low energy 
beta-beams |102j or spallation sources |103j . could of- 
fer a unique opportunity to precisely measure such cross 
sections and also realize fundamental interaction stud- 
ies, like measurements of the Weinberg angle at low mo- 
mentum transfer, a test of the Conserved- Vector-Current 
hypothesis (see [26] for a review), or of the existence of 
sterile neutrinos |31j . Note that improving the knowledge 
of the isospin and spin-isospin nuclear response involved 
in neutrino-nucleus interactions is of importance also for 
searches on the neutrino nature, since this constitutes a 
key ingredient of neutrino-less double-beta decay predic- 
tions [TU4] . 

Besides the fundamental particle or interaction proper- 
ties it encodes, observing a supernova luminosity curve is 
also of great astrophysical value. In fact, the time signal 
follows closely the different phases of the explosion, from 
the collapse to the accretion phase and to the cooling of 
the remaining proto-neutron star. This would provide 
key information to the understanding of the explosion 
mechanism of iron-core supernovae, a longstanding open 
problem. Currently, 2- and 3-dimensional simulations 
are being developed that realistically include the neu- 
trino transport, convection and hydrodynamical instabil- 
ities (the Standing- Accretion-Shock o SASI mode). Sev- 



eral groups are obtaining successful explosions for differ- 
ent progenitor masses, while reaching a consensus on the 
mechanism requires further studies (see e.g. [75l I105j ). 
Using 3D models based on a simplified neutrino trans- 
port scheme, Ref . [106] shows for example that one could 
be able to track the SASI mode in IceCube for a super- 
nova within 2 kpc. These massive stars are also candi- 
date sites for the r-process, whose identification is still an 
open issue; while several sites might contribute, such as 
the accretion disks around black holes [IB) . The impact 
of neutrino properties on the r-process is the object of 
numerous studies, as e.g. Ref. [107] that points out the 
possible role of sterile neutrinos in getting a successful 
nucleosynthesis in such sites. An intriguing question is 
also the possible impact of flavour conversion on the ex- 
plosion itself. MSW effects occur in the outer star layers 
and have no impact on the shock while collective effects 
from the neutrino-neutrino interaction take place deep in 
the star. The first investigations seem to indicate that 
they are still far out to affect the shock waves |108j . al- 
though it might still be too early to draw conclusions. 

A new window on the universe is opened by high- 
energy neutrino telescopes like ANTARES, the first un- 
dersea in the Mediterranean [109] . and IceCube, buried 
in deep ice at the South Pole |110) . The main mis- 
sion is to search for galactic and extra-galactic sources 
of high-energy neutrinos to elucidate the source of cos- 
mic rays and the astrophysical mechanisms that produce 
them, while these telescopes also investigate neutrino os- 
cillations, dark matter and supernova neutrinos (for Ice- 
Cube). At present, the IceCube collected data on ultra- 
high-energy neutrinos (above 10 19 eV) from active galac- 
tic nuclei (AGN) and gamma-ray-bursts (GRB) show no 
significant deviation from background in the number of 
observed events, although uncertainties in the models are 
still quite large. Interestingly two events have been re- 
cently identified at lower energy with a significance of 
about 3cr of not being a background |110j . High-energy 
neutrino telescopes are currently also providing data 
on neutrino oscillations measuring atmospheric neutri- 
nos, commonly a background for astrophysical neutrino 
searches. Using low energy samples, both ANTARES 
[111] and IceCube/DeepCore |112j have measured the pa- 
rameters 623 and Am23 in good agreement with existing 
data. Crucial information can be extracted from such 
telescopes in the future from low energy (below 40 GeV) 
atmospheric neutrinos. From the measurement of flavour 
conversion effects in the Earth predictions show that Ice- 
Cube/DeepCore offers a precise measurement of oscilla- 
tion parameters (see e.g. [1131 1114j ) and that PINGU, 
IceCube extension in the 10 GeV energy range, could 
measure the mass hierarchy and be sensitive to the Dirac 
phase |115j . Feasibility studies are currently ongoing 
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both for PINGU and for ORCA to establish if the en- 
ergy and angular resolution required for the mass hier- 
archy search can be achieved. Information on the phase 
could also be extracted by a precise measurement of high- 
energy neutrino flux ratios [116H118] . Ref . [118] gets a 2a 
coverage of about 10 % of the CP values, if uncertainties 
on the flux ratios are kept below 7 %. As for the fourth 
sterile neutrino its existence might lead to a distortion of 
the zenith angle distribution of high-energy atmospheric 
neutrinos through the MSW active-sterile resonance in- 
side the Earth. An analysis of data from Amanda and 
(incomplete) IceCube as well as the prospects for the 
completed IceCube give limits that for some of the sterile 
oscillation parameters are stronger than all combined ex- 
periments [119] . Ref. |120j shows that low energy events 
in DeepCore can substantially constrain the mixing of 
sterile neutrinos in the eV mass scale. The effect found 
is different for normal or inverted hierarchy of active neu- 
trinos, so that if such v s exist the hierarchy can also be 



identified. Finally both terrestrial experiments and neu- 
trino telescopes |121j set tight limits on Lorentz and CPT 
violation. A constraint also comes from the SN1987A 
from the nearly simultaneous arrival of the photons and 
neutrinos [122J. Interestingly, models based on Lorentz 
and CPT violation are developed to interpret oscillation 
data from all experiments (sec [123 and the nice sum- 
mary in Ref. [124] 1. 

In conclusion, neutrinos are messengers having a wide 
energy range, capable of covering cosmological distances 
and telling us about quiet and violent phenomena. Neu- 
trino physics is a domain rich of interdisciplinary as- 
pects and approaches. The measurement of these elu- 
sive particles requires inventive detectors of huge sizes. 
Sometimes slow on a man-scale, the progress continues 
steadily, while the last decade has been rich of excit- 
ing observations. Neutrinos have brought milestones in 
our understanding of fundamental issues in high-energy 
physics and astrophysics, and will likely bring more in 
the future. 



R. Davis, Phys. Rev. Lett. 12, 303 (1964). 

A. Suzuki, J. Phys. Conf. Ser. 120, 072001 (2008). 

K. Hirata et AL., Phys. Rev. Lett. 58, 1490 (1987). 

R. M. Bionta et AL., Phys. Rev. Lett. 58, 1494 

(1987). 

Alexeyev, E.N. et AL. , Proc. of the Leptonic Session 

of the 22nd Rencontre de Moriond (1987) 739. 

M. Aglietta et al. , Europhys. Lett. 3, 1315 (1987). 

G. Pagliaroli, F. Vissani, M. L. Costantini, and 
A. Ianni, Astropart. Phys. 31, 163 (2009). 

Y. Fukuda et AL., Phys. Rev. Lett. 81, 1562 (1998). 

H. A. Bethe, Phys. Rev.55, 1562 (1939). 

G. Gamow, and M. Schoenberg, Phys. Rev. 58, 1117 
(1940). 

B. Pontecorvo,, Sov. Phys. JETP 6, 429 (1957) [Zh. 
Eksp. Teor. Fiz.33, 549 (1957)]. 

Z. Mart, M. Nakagawa, and S. Sakata, Prog. Theor. 
Phys.28, 870 (1962). 



J. Beringer et al., Phys. 
C. Giunti, and C. W. Kim 
Physics and Astrophysics," 
p. 710. 

K. Abe et al., Phys. Rev. 
Y. Abe et al., Phys. Rev. 



Rev. D 86, 010001 (2012). 
"Fundamentals of Neutrino 
(Univ. Pr., Oxford, 2007), 



Lett. 107, 
Lett. 108. 



041801 (2011). 
131801 (2012). 



K. Ahn et AL., |arXiv:l 204.0626 [hep-ex]. 
P. An et al., Phy s. Rev. Let t. 108, 171803 (2012). 
Osipowicz et al., hep-ex/0109033 
Giunti and M. Laveder, Phys. Rev. C textbf83, 
065504 (2011). 

C. Athanassopoulos et al., Phys. Rev. Lett. 77, 
3082 (1996). 

A. A. Aguilar-Arevalo et al., Phys. Rev. Lett. 102, 



[23] 

[24' 
[25 : 
[26' 
[27 

[28' 

[29; 

[30 
[31 

[32' 
[33 

[34 

[35' 
[36 



LZ 



[37 

[38; 

[39 
[40 

[41 

[42' 
[43 



(2009) 101802. 

A. A. Aguilar-Arevalo et al., arXiv: 1207.4809 
[hep-ex 

G. Mention et al., Phys. Rev. D 83, 073006 (2011). 
P. Zucchelli, Phys. Lett. B 532, 166 (2002). 
C. Volpe, J. Phys. G 34, Rl (2007). 
A. Bandyopadhyay et al., Rept. Prog. Phys. 72, 
106201 (2009). 

J. Alonso et AL., jarXiv:1006.0260 
P. Huber, M. Lindner, T. Schwetz, and W. Win- 
ter, JHEP 0911, 044 (2009). 

M. Cribier et al., Phys. Rev. Lett. 83, 201801 (2011). 
C. Espinoza, R. Lazauskas, and C. Volpe, Phys. 
Rev. D 1986, , 113016 (2012). 



K. N. Abazajian et al. , |arXiv: 1204. 5379 [hep-ph]. 
J. LESGOURGUES and S. PASTOR, Adv. High Energy 
Phys. 2012, 608515 (2012). 

K. N. Abazajian et al., Astropart. Phys. 35, 177 
(2011). 

D. J. Koskinen, Mod. Phys. Lett. A 26, 2899 (2011). 



see http: / /agenda. infn.it /conferenceDisplay.py? 
confld=55 10| ~~ 
see http:/ /snews. bnl.gov/ 

J. F. Beacom, and M. R. Vagins, Phys. Rev. Lett. 93, 
171101 (2004). 

D. Autiero et al., JCAP 0711, 011 (2007). 

S. Calais, J. Kneller, C. Volpe, and J. Gava, Phys. 

Rev. D 81, 053002 (2010). 

J. F. Beacom, Ann. Rev. Nucl. Part. Sci. 60, 439 
(2010). 

C. Lunardini, |arXiv:1007.3252| [astro-ph.CO], 
L. Wolfenstein, Phys. Rev. D 17, 2369 (1978). 



10 



S. P. Mikheev, and A. Y. Smirnov, Sov. J. Nucl. [76 
Phys.42, 913 (1985). 

Q. R. Ahmad et al., Phys. Rev. Lett. 89, 011301 [77 
(2002). 

K. Eguchi et al., Phys. Rev. Lett. 90, 021802 (2003). [78 

G. Bellini et al., Phys. Rev. Lett. 108, 051302 
(2012). [79 
A. M. Dziewonski, and D. L. Anderson, Phys. Earth 
Planet. Interiors 25, 297 (1981). [80 
V. K. Ermilova, V. A. Tsarev, and V. A. Chechin, [81 
JETP Lett. 43, 453 (1986) [Pisma Zh. Eksp. Teor. 
Fiz.43, 353 (1986)]. [82 
E. Kh. Akhmedov, Yad. Fiz. 47, (1988) 475 [Sov. J. [83 
Nucl. Phys. 47 (1988) 301]. [i 
E. Kh. Akhmedov, Nucl. Phys. B 538, 25 (1999). [85 
S. T. Petcov, Phys. Lett. B 434, 321 (1998). [86^ 
E. K. Akhmedov, A. Dighe, P. Lipari, and 

A. Y. Smirnov, Nucl. Phys. B 542, , 3 (1999). [87 
E. K. Akhmedov, M. Maltoni, and A. Y. Smirnov, 
JHEP 0806, 072 (2008). [. 

A. S. Dighe, and A. Y. Smirnov, Phys. Rev. D 62, 

033007 (2000). [. 

R. C. Schirato, and G. M. Fuller, |aitro- ph/0205390"l 

[90 

R. Tomas et al., JCAP 0409, 015 (2004). 

B. Dasgupta, and A. Dighe, Phys. Rev. D 75, 093002 [91 

(2007) . 

F. N. Loreti et AL., Phys. Rev. D 52, 6664 (1995). [92 

G. L. Fogli, E. Lisi, A. Mirizzi, and D. Mon- 

tanino, JCAP 0606, 012 (2006). [93 

A. Friedland, and A. Gruzinov, |astro-ph /0607244 
J. P. Kneller, and C. Volpe, Phys. Rev. D 82, [94 
123004 (2010). [95 

H. Duan, and J. P. Kneller, J. Phys. G 36, , 113201 
(2009). [96 

J. T. Pantaleone, Phys. Lett. B287, 128 (1992). 

S. Samuel, Phys. Rev. D 48, 1462 (1993). 
A. B. Balantekin, and H. Yuksel, New J. Phys. 7, [97 
51 (2005). 

H. Duan, G. M. Fuller, J. Carlson, Y.-Z. Qian, [i 
Phys. Rev. D 74, 105014 (2006). [99 
H. Duan, and G. M. Fuller, and Y. -Z. Qian, Ann. 
Rev. Nucl. Part. Sci.60, 569 (2010). [100 

H. Duan, G. M. Fuller, J. Carlson, and Y. - [101 
Z. Qian, Phys. Rev. D 75, .125005 (2007). 

S. Hannestad, G. G. Raffelt, G. Sigl , and [102 
Y. Y. Y. WONG, Phys. Rev. D 74, 105010 (2006) [103 
[Erratum-ibid. Phys. Rev. D 76, 029901 (2007)]. 

G. G. Raffelt, and A. Y. Smirnov, Phys. Rev. D 76, [104 
081301 (2007) [Erratum-ibid. Phys. Rev. D 77, 029903 [105 

(2008) ]. [106 : 
S. Calais, and C. Volpe, Phys. Rev. D 84, 085005 [107 
(2011). 

A. Esteban-Pretel et al., Phys. Rev. D 78, 085012 
(2008). [108 
S. Chakraborty et al., Phys. Rev. Lett. 107, 151101 
(2011). [109 

H. -T. Janka, Ann. Rev. Nucl. Part. Sci.62, 407 (2012). 



and C. Espinoza. 



J. F. Cherry et al., Phys. Rev. Lett. 108, 261104 
(2012). 

A. B. Balantekin, and Y. Pehlivan, J. Phys. G 34, 
47 (2007). 

C. Volpe, D. Vaananen, 
larXiv:1302.2374l 

R. Surman, G. C. McLaughlin, and W. R. Hix, As- 
trophys. J. 643, 1057 (2006). 

L. -T. Kizivat et AL., Phys. Rev. C 81, 25802 (2010). 
A. Malkus, J. P. Kneller, . C. McLaughlin, and 
R. Surman, Phys. Rev. D 86, 085015 (2012). 
A. D. DOLGOV, Phys. Rept.370, 333 (2002). 

F. IOCCO et AL., Phys. Rept.472, 1 (2009). 

J. Hamann et al., Phys. Rev. Lett. 105, 181301 (2010). 
S. Weinberg Phys. Rev. 128, 1457 (1962). 
A. G. Cocco, G. Mangano, and M. Messina, 
JCAP0706, 015 (2007). 

R. Lazauskas, P. Vogel, and C. Volpe, J. Phys. G 
35, 025001 (2008). 

S. Davidson, E. Nardi, and Y. Nir, Phys. Rept.466, 
105 (2008). 

H. Minakata, and S. Watanabe, Phys. Lett. B468, 
256 (1999). 

E. K. Akhmedov, C. Lunardini, and A. Y. Smirnov, 
Nucl. Phys. B 643, 339 (2002). 

A. B. Balantekin, J. Gava, and C. Volpe, Phys. 
Lett. B662, 396 (2008). 

J. P. Kneller, and G. C. McLaughlin, Phys. Rev. 
D80, 053002 (2009). 

J. Cava, and C. Volpe, Phys. Rev. D 78, 083007 

(2008) . 

J. Cava, and C. Volpe, Nucl. Phys. B 837, 50 (2010). 
V. Barger, P. Huber, and D. Marfatia, Phys. Lett. 
B 617, 167 (2005). 

J. Cava, J. Kneller, C. Volpe, and 

G. C. McLaughlin Phys. Rev. Lett. 103, 071101 

(2009) . 

A. S. Dighe, M. T. Keil, and G. G. Raffelt, JCAP 
0306, 005 (2003). 

P. D. Serpico et al., Phys. Rev. D 85, 085031 (2012) 

D. Vaananen, and C. Volpe, JCAP 1110, 019 

(2011) . 

ht t p : / / www . phy. duke.edu/~schol/ snowglobes / 

K. Scholberg, Ann. Rev. Nucl. Part. Sci.62, 81 

(2012) . 

C. Volpe, J. Phys. G 30, LI (2004). 

F. T. Avignone, and Y. .V. Efremenko, J. Phys. G 
29, 2615 (2003). 

C. Volpe, J. Phys. G 31, 903 (2005). 

A. Burrows, Rev. Mod. Phys. 85, 245 (2013). 

T. Lund et al., Phys. Rev. D 86, 105031 (2012). 

G. C. McLaughlin, J. M. Fetter, A. B. Bal- 
antekin, and G. M. Fuller, Phys. Rev. C 59, 2873 
(1999). 

B. Dasgupta, E. P. O' Connor and C. D . Ott, Phys. Rev. 
D 85, 065008 (2012) [arXiv: 1106. 1167] [astro-ph. SR] ] . 
M. Ageron et al., Nucl. Instrum. Meth. A 656, 11 
(2011). 



11 



[110] A. Kappes, Phys. Conf. Ser. 09, 012014 (2013). 
[Ill] S. Adrian-Martinez et al., Phys. Lett. B714, 224 
(2012). 

[112] A. Gross, |arXiv:1301.4339| [hep-ex], 

[113] E. Fernandez-Martinez, G. Giordano, O. Mena, 

and I. MOCIOIU, Phys. Rev. D 82, 093011 (2010). 
[114] O. Mena, I. Mocioiu, and S. Razzaque, Phys. Rev. 

D 78, 093003 (2008). 
[115] E. K. Akhmedov, S. Razzaque, and A. Y. Smirnov, 

JHEP 02, 082 (2013) [JHEP 1302, 082 (2013)] 
[116] P. D. Serpico, and M. Kachelriess, Phys. Rev. Lett. 

94, 211102 (2005). 



[117] W. Winter, Phys. Rev. D 74, 033015 (2006). 
[118] D. Meloni, and T. Ohlsson , Phys. Rev. D 86, 
067701 (2012). 

[119] A. Esmaili, F. Halzen, and O. L. G. Peres, JCAP 

1211, 041 (2012). 
[120] S. Razzaque, and A. Y. Smirnov, Phys. Rev. D 85, 

093010 (2012). 

[121] R. Abbasi et al., Phys. Rev. D 82, 112003 (2010). 
[122] M. J. LONG O, Phys. Rev. D 3 6, 3276 (1987). 
[123] T. Kator i, |arXiv:1211.7129| [hep-ph], 
[124] J. S. Diaz, |arXiv: 1109.4 620 [hep-ph]. 



